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1 Introduction 

We consider one dimensional solutions and higher dimensional radial solutions of 
the following Klein-Gordon equation: 



d 2 t U = AU + lUl^U - U, 
U(0) = U and U t (0) = U u 

where U(t) : x G R N -> U(x,t) G E, U Q G and XJ X G l£ C)U . 

The space Lf oc u is the set of all v in Lf oc such that 

\\v\\ L 2 = sup I / |f(2;)| 2 da; I < +oo, 

loC ' U aeR" \J\x-a\<l J 

and the space B.} oc u = {v \ v, Vv G Lf oc u }. 



'i.r 



The nonlinear Klein Gordon equation appears as a model of self-focusing waves 
in nonlinear optics (see Bizon, Chamj and Szpak [4]). 



More generally, we consider the following semilinear wave equation: 

d 2 U = AU+ \ur l U + f(U) + g(\x\,t, Vt/.g, d t U), 
U(0) = U and U t (0) = U x . 



;i.2) 



We assume that the functions / and g are C l functions, where / : R — > R and 
g : R 4 — > R satisfy the following conditions 

(H f ) \f{u)\ < M(l + \u\ q ), for all y eR with (q < p, M > 0), 

(H g ) \g(X,t,v,z)\< M(l + \v\ + \z\), for all X, t, v, z e R with (M > 0). 

We assume also that the function g is globally Lipschitz. Finally, we assume that 

4 

l<pandn<l + if JV > 2. (1.3) 

N - 1 - v 1 

Since J7 is radial if > 2, we introduce 

u(r, t) = U(r, t) for r e R, if JV = 1 (1.4) 

and 

u(r, t) = U(x,t) if r = |x| and > 2 (1.5) 

and rewrite (1.2) as 

d\u = d^u + ^ N ~^ d r u + |n| p_1 -u + f(u) + g(r, t, d r u, d t u), 
d r u(0,t) = if N > 2, (1.6) 
w(r, 0) = uq{t) and u t (r, 0) = Ui(r), 

where w(t) : r G / -»■ u(r, t) G R, with / = R+ if N > 2 and / = R if N = 1. 

The Cauchy problem of equation (1.2) is solved in H} ocu x Lf ocu . This follows 
from the finite speed of propagation and the wellposedness in H 1 x L 2 (see for 
example Georgiev and Todorova [8]) , valid whenever 1 < p < 1 + jr^o- The 
existence of blow-up solutions u(t) of (1.2) follows from energy techniques (see for 
example Levine and Todorova [15]) and Todorova [27]). 

If u is a blow-up solution of (1.6), we define (see for example Alinhac [1]) a 1- 
Lipschitz curve T = {(r, T(r))} where r G / such that the maximal influence domain 
D of u (or the domain of definition of it) is written as 

D = {(r,t) \ re I, t< T(r)}. (1.7) 



2 



r is called the blow-up graph of u. A point tq > is a non-characteristic point if 
there are 

5 G (0, 1) and t < T(r ) such that u is defined on C ro T i ro \ $ H {t > t } fl {r G /} 

(1.8) 

where C^f j = {(r,t) \ t < t — S\r — r|}. We denote by 1Z C / (resp. S C I) the set 
of non-characteristic (resp. characteristic) points. 

In the case (/, g) = (0, 0), equation (1.2) reduces to the semilinear wave equation: 

r %u = au + lur'u, 

\ U(0) = U and U t (0) = U x . K ' 

In a series of papers [20], [21], [24] and [22] (see also the note [23]), Merle and 
Zaag give a full picture of the blow-up for solutions of (1.9) in one space dimension. 
Recently, in [7], Cote and Zaag refine some of those results and construct a blow- 
up solution with a characteristic point a, such that the asymptotic behavior of the 
solution near (a,T(a)) shows a decoupled sum of k solitons with alternate signs. 
Moreover, in [25], Merle and Zaag extend all their results to higher dimensions in 
the radial case, outside the origin. Our aim in this work is to generalize the result 
obtained for equation (1.9) in [25] to equation (1.2). Let us note that all our results 
and proofs hold in both cases of (1.6) (N = 1 and N > 2). However, the situation 
is a bit more delicate when N >2, since we have to avoid the origin which brings a 
singular term ^^-d r u to (1.6). Thus, for completeness, we focus on the case N > 2 
avoiding r = 0, and stress the fact that all our results hold in the case N — 1, even 
when r = 0, and with no symmetry assumptions. 

Throughout this paper, we consider U(x, t) a radial blow-up solution of equation 
(1.2), and use the notation u(r,t) introduced in (1.4). We proceed in 3 sections: 

- in Section 2, we give a new Lyapunov functional for equation (1.6) and bound the 
solution in the energy space. 

- in Section 3, we study 71, in particular the blow-up behavior of the solution and 
the regularity of the blow-up set there. 

- in Section 4, we focus on S, from the point of view of the blow-up behavior, the 
regularity of the blow-up set and the construction of a multi-soliton solution. 

We are aware that our analysis is a generalization of the radial case of equation 
(1.9) treated by Merle and Zaag in [25]. For that reason, we will give the statements 
of the results for equation (1.2) and focus only on how to deal with the new pertur- 
bation terms appearing in (1.2). Let us add that, we believe that our contribution 
is non trivial and introduces a new approach for perturbed problems. Moreover, it 
proves a bunch of results, especially for the Klein-Gordon equation (1.1). 
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2 A Lyapunov functional for equation (2.2) and a 
blow-up criterion in the radial case 



We showed in [11] and [10] that the argument of Antonini, Merle and Zaag in 
[3], [17], [19] and [18] extends through a perturbation method to equation (1.2) 
with no gradient terms even for non radial solutions. The key idea is to modify 
the Lyapunov functional of [3] with exponentially small terms and define a new 
functional which is decreasing in time and gives a blow-up criterion. In [25], Merle 
and Zaag successfully used our ideas to derive a Lyapunov functional for the radial 
case with no perturbations (i.e. for equation (1.6) with (/, g) = (0,0)). Here, we 
further refine our argument in [11] and [10] to derive a Lyapunov functional for 
equations bearing the two features: the presence of perturbation terms and the 
radial symmetry. For the reader's convenience, we briefly recall the argument in the 
following. 

Given ro > 0, we recall the following similarity variables' transformation 

w ro (y,s) = (T(r )-t)^u(r,t), V = Z~s°_ v s = -log(T(r ) - t). (2.1) 

The function w = w ro satisfies the following equation for all y G (—1,1) 
and s > max (— log T(r ), — logr ): 

P2 r 2 (P +1 ) ,| ip-l P + 3 x o ^2 

OeW = Lw — --w + \wr w o s w — 2yo w 

(p — I) 1 p — 1 y ' 

(N — 1) 2ps / 2s \ , 

+e~ s ^ ^d y w + e-—^f[e—^w) (2.2) 



r + ye~ 

2ps / _ _ (p + l)s (p+l)s 2 

+e p- 1 g r + ye s ,T — e s ,e p- 1 d v w,e p- 1 (d s w + yd v w H mi' 

V p — 1 ' 

1 2 

where Cw = -d y (p(l — y 2 )d y w) and p(y) = (1 — y 2 )p- 1 . (2.3) 
In the whole paper, we denote 

F{u) = [ f(v)dv. (2.4) 



Let us recall that for the case (/,<?) = (0,0), the Lyapunov functional in one 
space dimension is 

Eo(w(s)) = £ {^-{d s wf + \ (d y wf (1 - y 2 ) + j^yy ~ ^H^pdy, 

(2.5) 
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which is defined in the Hilbert space 

K=l [ qe HL x LL(-1, 1) I \\qfu = J ^ (ql + (q[f (1 - y 2 ) + ql) pdy < +00} . 

(2.6) 

Introducing 

E(w(s),s) = E (w(s)) + I(w(s),s) + J(w(s),s) (2.7) 

where, 

/1 
F(ep- 1 w)pdy, (2.8) 

J(w(s),s) = -e~ ls [ wd s wpdy, (2.9) 



1 



with 



we claim the following: 



7 = min(^ £—2) >0, (2.10) 
2 p — 1 



Proposition 2.1 (A new functional for equation (2.2)) 

(%) There exists C = C(p, N, M,q) > and Sq(p, N, M, q) £ R snc/i that for all 
r > and /or a// s > max (— logT(r ), So, — 41ogr , — log yV 

as z p 1 j _i 1 — y 

(2.11) 

(^zj (A blow-up criterion) There exists Si(p, N, M, q) £ R such that, for all 
s > max (s 0) £1, ), we /icrae H(w(s), s) > 0. 

Remark: From (i), we see that the Lyapunov functional for equation (2.2) is in 
fact H(w ro (s),s) where 

H(w(s),s) = E(w(s),s)e E ir e "' s + pe~^ s (2.12) 
not E(w ro (s), s) nor E (w ro (s)), for some large constant p. 

Remark: We already know from [11] and [10] that even in the non-radial setting, 
equation (1.6) has a Lyapunov functional given by a perturbed form of a natural 
extension to higher dimensions of E(w ro (s), s) (2.5). Unfortunately, already for the 
non-perturbed case of (1.2) with (f,g) = (0,0), due to the lack of information on 
stationary solutions in similarity variables in dimensions N > 2, it wasn't possible to 
go further in the analysis, and the investigation had to stop at the step of bounding 
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the solution in similarity variables. On the contrary, when N = 1, Merle and Zaag 
could obtain a very precise characterization of blow-up in [20], [21], [24], [22] (with 
some refinements by Cote and Zaag in [7]). 

Here, considering perturbations as stated in (1.2) and restricting ourselves to one 
dimensional solutions or higher dimensional radial solutions, we find a different 
Lyapunov functional. Considering arbitrary blow-up points in one-space dimen- 
sion(including the origin) and any non-zero blow-up in higher dimensions, the char- 
acterization of stationary solutions in one space dimension is enough, and we are 
able to go in our analysis as far as in the one-dimensional case. 

Following [17] and [20], together with our techniques to handle perturbations in 
[11] and [10], we derive with no difficulty the following: 

Proposition 2.2 (Boundedness of the solutions of equation (2.2) in the 
energy space) For all r > 0, there is a C^ro) > and S^^o) £ M such that for 
all re and s > S 2 (r ), 

J' ((d y w r ) 2 (l - y 2 ) + K) 2 + {d s w r f + Kr 1 ) pdy < C 2 (r ). 

Proof. The adaptation is straightforward from [17] and Proposition 3.5 page 66 in 
[20]. The only difference is in the justification of the limit at infinity of Eo(w ro (s)) , 
which follows from the limit of H(w ro (s), s) defined in (2.12). In fact, we know from 
Proposition 2.1 that H(w ro (s),s) is decreasing and bounded from below, and such 
an information is unavailable for E (w ro (s)). ■ 

Proof of Proposition 2. 1: 
(i) Consider r > 0, s > max(— log T(r ), 0, — log y, — 41ogr ) and write w = w ro 
for simplicity. From the similarity variables' transformation (2.1), we see that 



r + ye £ 



ro 3ro 
2 ' 2 



(2.13) 



Multiplying equation (2.2) by d s wp and integrating for y £ (—1,1), we see by (2.5) 
and (2.8) that 

±(E (w(s)) + I(w(s),s)) = ^- T i^l p dy+(N-l)e- s f d s wd y w% 
ds x p - 1 J_ x 1 - y 2 J_! r 



h(s) 

2(p+l) 2(p+l)s f 1 / _2s_ \ 2 _2p± f 1 ( _2s_ \ 

H e p- 1 / F[ep- 1 w)pdy -\ e p- 1 / f \ep- 1 w)wp<ly 

P-l J -i v ' P~ l J-i ^ ' 

V v ' V v ' 

h(s) h(s) 

2ps f 1 / _ _ (p + l)s (p+l)s 2lU ,\ 

+e p- 1 / g[r Q + ye s ,T — e s ,e p- 1 d y w,e p- 1 (d s w + yd y w H -) j d s wpdy . 

J -i ^ p — 1 / 



h(s) 
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where r is defined in (2.13). Using (2.13), we write 

\h(s)\ <Ce~ s [\d y w) 2 p(l-y 2 )dy + ^- [ (d s wf— P —dy. (2.14) 

Using the fact that 

\F{x)\ + |x/(x)| < C(l + \x\ q+1 ) < C(l + |x| p+1 ), (2.15) 

where F and / are defined in (2.4) and (1.2), we obtain that 

»i 

|/ 2 (s)| + \I 3 (s)\ < Ce~ 2J l^ + Ce- 2J l^ J \w\ p+1 pdy. (2.16) 
Using the inequality ab < y- + y and the hypothesis (i/ 9 ), we write that 

\h{s)\<Ce~ s j ({d s w) 2 + (w 2 )^)pdy + Ce- s j \d y w\\d s w\pdy + Ce~ s . (2.17) 

Similarly, we prove that 

^ I^H^Ipd?/ < J 1 {d.wfj^Ay + j\d y wf{l - y 2 )pdy. (2.18) 

Combining (2.17) and (2.18), we conclude that 

I < Ce~ s ^({dywfil - \y\ 2 ) + ^£ + W 2 )pdy + Ce~ s . (2.19) 

Then, by using (??), (2.14), (2.16) and (2.19), we deduce that 

±(E (w(s)) + I(w(s) lS )) < ( i- + Ce-i) (d s wf-^dy 

as p — 1 7-i 1 — y 

+Ce -2~fs j \ w \p+i pdy + Ce- 2 ^ s . (2.20) 

Considering J(w(s),s) defined in (2.9), we obtain from equation (2.2) and integra- 
tion by parts 

els is~ JMs) ' s) = ~ / (MVi/ + Jj 9 ^ 1 - y 2 )p^ + t^tt^ / w2 p d y 

/ 1 p _|_ 3 r l r 1 

\w\ p+1 pdy + (7 H 2AH / wd s w(s)pdy — 2 / wd s wyp'dy 
-1 P-l 7-i 7-1 

" ^ J d s wdywypdy — e~~ J wf(e~w S jpdy — (N—l)e~ s J wd y w-dy 

2ps f 1 / __ _ (p + l)s (p+l)s 2 \ 

— e p_1 y ^(jo + ye S ,T — e s ,e p- 1 d y w,e p- 1 (d s w + yd y w-\ -wjjpdy. 



Combining (2.7), (2.8) and (??), we write 

ev£j(w(s),s)< p -f (E (w(s)) + I(w(s), a))-*? / {d y wf{\ - y 2 )pdy 



1 

1 



+ (7 + ^| - 2N + ^-e^ s ) I wd s wpdy (2.21) 



Ms) 

8 f 1 y 2 f 1 _2e± f 1 ( _2*_ \ 

H / ludgtu -pdy—2 / d s wd v wypdy —e p- 1 / w/( ep- 1 ^ Ipdy 

p-lj-i i-y 2 J-i J-x \ / 

S v " v " v ' 

J 2 (s) J 3 (s) J4(s) 

2ps f 1 f - — + (p+l)s 2 \ 

; p l J w 9yo + y e s ,To — e s ,e p- 1 d y w,e p- 1 (d s w + yd y w H -wjjpdy 



V + 3 2(p+l)s Z" 1 2 n 

+ - e p- 1 / F(ep- lS w)pdy-(N -l)e~ s wd y w~dy . 

2 7-i J-i r 



v " v 

Je(s) J7{s) 



We now study each of the last five terms. To estimate Ji(s) and his), we use the 
Cauchy-Schwartz inequality to have 



< Ce 1 ^ J 1 id s w) 2 ^-^dy + Ce-^ J 1 w 2 pdy. (2.22) 
\Us)\ < C ^ f\dsw) 2 — P —dy + Ce-^ f w *-^dy. 

J-i 1 - y J-i 1 - y 

Recalling the following Hardy-Sobolev estimate (see Appendix B page 1163 in [17] 
for the proof): 

f h 2 -?— 2 dy <C f h 2 pdy + C f ih'iy)) 2 pil - y 2 )dy, (2.23) 

J-i 1 — y J-i J -i 



we conclude that 



\Us)\ < Ce^ J' id s w) 2 y^- 2 dy + Ce-^ J 1 w 2 pdy 

+Ce-^ J\d y w) 2 pil - y 2 )dy. (2.24) 
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Using the Cauchy-Schwartz inequality, we have 

I JsOOl < Ce? Jjdswf^^&y + Ce~^ j\d y wfp{\ - y 2 )dy. 
From (2.15), we write 

|J 4 (s)| + |J 6 (s)| < Ce^ s + Ce^ s J' \w\ p+1 pdy. 
In a similar way, using the hypothesis (H g ) and (2.23), we have 



(2.25) 



(2.26) 



-7« / fan„\z (L 

1 - y* 



\Ms)\ < Ce^ s I (dswY^—^y + Ce-^ (d y w) 2 p(l - y 2 )dy 



1 

2 



+Ce~~< s / w 2 pdy + Ce^ s . (2.27) 



1 



Using (2.13) and (2.23), we write 

\J 7 {s)\<Ce~i J {d y w) 2 p{\ - y 2 )dy + Ce'i J w 2 pdy. (2.28) 
Finally, using (2.21), (2.22), (2.24), (2.25), (2.27) and (2.28), we deduce that 

e^-^J(w(s),s) < P±^(e q (w(s))+I(w(s),s)) (2.29) 

+ {Ce^ - ^) J\d y w) 2 (l - y 2 )pdy 



p+1 



1 2 p 



+Ce^ / (d s w) 2 -^dy + Ce^ s . 
1 - y 2 



From (2.20) and (2.29), we obtain 

. P + 3 

[w{s ,s 



d -E(w(s),s) < Ce'*" + E±^ e --r»E(i 



ds v v " ' ~ 2 



Ce-^--^) f {d s w) 2 -^- 2 dy 



p-l/ J-! l-y z 



2(p-l), 

+ (Ce- : ¥ - P -^)e^ s f\d y wf{l - \y\ 2 )pdy (2.30) 
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We now choose So > 0, large enough, so that for all s > So, we have 



P^l_Ce-^>0, -l±l--Ce~ : f>0, -l^L _ Ce~^ > 0, -^—-Ce^X). 
4 2(p-l) - ' 2(p + l) ~ p-1 

Then, we deduce that, for all s > max(So, — logT(r ), — log 2 ^, — 41ogr ), we have 



£e{w(s),8) < Ce-*" + P±^ e -r^(«,( a ), s) - ^ j\d s w) 2 ^^dy . (2.31] 



This yields (i) of Proposition 2.1. 

(ii) We finish the proof of Proposition 2.1 here. More precisely, we prove that 
there exists S\(p, N, M, q) G M such that, for all x £ and T G (0, T(x )], 

Vs> max (-log T , S u ) H(w xo , To (s), s) > 0. (2.32) 

We give the proof only in the case where Xq is a non characteristic point. Note 
that the case where Xq is a characteristic point can be done exactly as in Appendix 
A page 119 in [20]. If x is a non characteristic point, the argument is the same 
as in the corresponding part in [3]. We write the proof for completeness. Arguing 
by contradiction, we assume that there exists a non characteristic point xq G Mr, 
T G (0,T(x o )] and s\ > max (—log To, Si, ) such that H(w(si), s±) < 0, where 
w = w Xq ^t - By definition (2.12) of H, we write 

H(W(s), s) > /le" 27 * - e¥ e ~~' a (— ^- + Ce~ 2 ^ J \W\ p+1 pdy 



l / „ , 1 \ rl 



+e£ +3 e _ 7S / / 1 Ce \ r {dsW) 2 pdy + ( P + i _ Ce ^ s j i W 2 pdy 



1 



P+ 1 



> / \W\ p+1 pdy. 



i 



-i 



if s > S 2 (p,N,q,M) > S!(p,N,q,M) for some S 2 (p,N,q,M) G M large enough. 
Using this inequality together with the fact that H(W(s),s) is decreasing by the 
remark following Proposition 2.1, we see that the argument used by Antonini and 
Merle in Theorem 2 page 1147 in [3] for the equation (1.9) works here and we get 
the blow-up criterion. This concludes the proof of Proposition 2.1. ■ 

3 Blow-up results related to non-characteristic points 

Let us first introduce for all \d\ < 1 the following solitons defined by 

(1 _ d?)*^ 1 f2(p+l)\p^ T 
K(d,y) = k,q — where k = I -, -r^ I and \y\ < 1. (3.1) 



dy)^ KiP- 1 )' 
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Note that n(d) is a stationary solution of (2.2), in the particular case where (/, g) = 
(0,0) and in one space dimension. 

Adapting the analysis of [20] and [21], we claim the following: 

Theorem 1 (Blow-up behavior and regularity of the blow-up set on TV) 

(i) (Regularity related to TV) 1Z ^ 0, 7£nlR+ is an open set, and x h> T(x) is of 
class C 1 on1ZnR* + . 

(ii) (Blow-up behavior in similarity variables) There exist fi > and Cq > 

such that for all r £ 1Z PI there exist 0(r ) = ±1 and s (r ) > — log T(r ) such 
that for all s > So-' 

|(ait))-«(o (T ' M) ) 

Moreover, E (w ro (s)) — > E (k ) as s — > oo. 

Remark: As stated in the introduction, this result holds also when N = 1, with no 
symmetry assumptions an initial data, for all ro £ M, even ro = 0; when N > 2 and 
if £ 1Z, the asymptotic behavior of remains open. 

Proof. As in the non-perturbed radial case (take (/, g) = (0,0) in (1.6)) treated in 
[25], we need to make some minor adaptations to the one-dimensional non-perturbed 
case treated in [20] and [21]. It happens that the same adaptation pattern works 
in the present case, and that is the reason why we don't mention it, and refer the 
reader to the proof of Theorem 1 in page 358 of [25]. The only points to check are 
the following: 

- Continuity with respect to the scaling parameter. Due to the fact that equation 
(1.6) is no longer invariant under the scaling 

A h-> u x (£, t) = A^w(A{, Ar), 

we need to understand the continuous dependence of the solutions of the following 
family of equations 

d 2 t u = dlu + ( N ~y X d r ii + lul^u + X^f( X^u) 
1 r x + Xr 11 V / 

+A^#(V, At, X~^d r u, X~^d t uj , (3.3) 

with respect to initial data and the parameters x > and A > (including the limit 
as A —7- and this is a classical estimate. 

- A new statement for the trapping result: This is due to the fact that equation 
(2.2) in similarity variables depends on a parameter ro > and contains new terms 
of order e~ ls (7 is defined in (2.10)) (it is no longer autonomous). This is the 
trapping result in our setting: 



< C e- w{s - so) . (3.2) 
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Theorem 2 (Trapping near the set of non zero stationary solutions of 

(2.2)) 

For all po > 0, there exist positive €q, po and Cq such that for all e* < €q, there 
exists s (e*) such that if r > p , s* > s and w G C([s*,oo),?i) is a solution of 
equation (2.2) with 



Vs > s* , E(w(s), s) > Eq(kq) — e 



7 s 



(3.4) 



and 



w{s*) 
d*w(s* 



to 







< e* 



H 



for some d* G (—1,1) and u* = ±1, then there exists d^ G (—1,1) such that 

largtanhtioo — argtanhcf | < Cq€* , 

and for all s > s* , 



w(s) 
d„w(s) 



— CO 



^idoo-, ■) 




H 



Proof. The proof follows the pattern of the radial case treated in [25]. For that 
reason, we refer the reader to the Proof of Theorem 2 page 360 in that paper, 
and focus in the following only on how to treat the new terms generated by the 
perturbations / and g in (1.2). With respect to the pure power case in one space 
dimension, the difference comes from the linearization of (2.2) around the stationary 
solutions K,(d,y) in (3.1), where we see the following lower order terms: 



(N — l)e" 



r + ye~ 



-dyW 



2ps 

e p- 1 

_ 2ps 

e p- 1 



f ( ep- x w 



< — (N — l)e~ s \d y w\; 

2(p-q)s 2(p-q)s 



(3.5) 



w\ 



_ _ (P+I)s (p+l)s 

g[r + ye s ,T -e s ,e p- 1 d y w,e p- 1 (d s w + yd y w 



-w 



< CMe 



1 + \d s w\ + | d y w | + \wn. 



P 



as soon as ro > po > 0, and s > — log^ 1 - F° r more details an the adaptation, we 
refer the reader to the proof of Theorem 1 in page 358 of 



4 Blow-up results related to characteristic points 

The first question in this case is of course the existence of examples of initial data 
with S ^ 0. If the perturbation g introduced in (1.2) does not depend on |x|, then 
the existence of such an example follows from the knowledge of the blow-up behav- 
ior at non-characteristic points, as in the pure power nonlinearity case (1.9). If g 



12 



depends on then we need to apply the constructive method of Cote and Zaag 
[7], which relies fundamentally on the knowledge of the blow-up behavior near a 
characteristic point. For that reason, we leave the existence issues to the end of the 
section, and start with the description of the blow-up features near characteristic 
points. More precisely, we proceed in two sections: 

- In Section 4.1, we consider arbitrary blow-up solutions having a non-zero charac- 
teristic point, and we give a full description of its blow-up behavior and its blow-up 
set near this characteristic point. 

- In Section 4.2, we prove the existence of such a solution, and also give some criteria 
for the existence or the non-existence of characteristic points. 

4.1 Description of the blow-up behavior and the blow-up 
set near a characteristic point 

Now, given r G S fl W + , we have the same description for the asymptotic of w ro 
as in the one- dimensional case with no perturbations (i.e. for equation (1.6) with 
(/, g) = (0, 0)) refined recently by Cote and Zaag in [7]. In order to state the result, 
let us introduce 

500= (i--^^)^-^ log* + (4.1) 

where the sequence (ai)i=i,...,fc is uniquely determined by the fact that (Ci( s ))i=i,...,fc 
is an explicit solution with zero center of mass for the following ODE system: 

_q _ g-jrr(Ci-Ci-i) _ g-^rrCCi+i-Ci^ ^2) 

Cl 

where c\ = c\(p) > and Co( s ) = Cfc+i( s ) = (see Section 2 in [7] for a proof of 
this fact). Note that c\ = c\{jp) > is a constant appearing in system (4.11), itself 
inherited from Proposition 3.2 of [24]. With this definition, we can state our result 
(for the statement in one space dimension, see Theorem 6 in [24] and Theorem 1 in 
[7]): 

Theorem 3 (Description of the behavior of w ro where r is characteristic) 

Consider tq e S fl Ml. Then, there is Co( r o) £ R such that 

— > and E (w ro (s)) —> k(r )E (K ) 

(4.3) 

fc(r ) > 2, (4.4) 
9 i ^e 1 (-l) i+ \ 0! = ±1 (4.5) 



w ro {s) 
d s w ro (s) 



0i 



i=i 



as s — > oo ; for some 
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and continuous di(s) 



tanh(j(s) with 



Ci(s) = Us) + Co 



(4.6) 



where Q(s) is introduced above in (4.1). 

Remark: As stated in the introduction, this result holds also when N — 1, with no 
symmetry assumptions an initial data, for all ro G M, even ro = 0; when N > 2 and 
if G S, the asymptotic behavior of wq remains open. 

Proof. As in the one-dimensional case with no perturbations (i.e. for equation (1.6) 
with (/, g) = (0, 0)), the proof of the asymptotic behavior and the geometric results 
on S (see Theorem 5 below) go side by side. Evidently the refined description given 
by (4.6) is obtained as in [7]. We leave the proof after the statement of Theorem 5. 

■ 

Let us note that we get the following result on the energy behavior from the 
asymptotic behavior at a non-characteristic point (see (ii) of Theorem 1) and at a 
characteristic point (see Theorem 3): 

Corollary 4 (A criterion for non-characteristic points) 

For all r > 0, there exist C3(r ) > and S3 (r ) G R such that: 
(i) For all r G [y, ^r] and s > S3, we have 



then r £ 1Z. 

Remark: With respect to the statement in one-space dimensions with no perturba- 
tions, (Corollary 7 in [24]), this statement has additional exponentially small terms. 
This comes from the fact that the functional E(w(s), s) is no longer decreasing, and 
that one has to work instead with the functional H(w(s), s) (2.12) which is decreas- 
ing, and differs from E(w(s),s) by exponentially small terms, uniformly controlled 
for r G [y, ^y-] thanks to the uniform estimates of Proposition 2.2. 
Proof. If one replaces E(w(s),s) by H(w(s), s), then the proof is straightforward 
from Theorems 1 and 3 together with the monotonicity of H(w(s), s) (see (2.12) and 
(2.1)). Since the difference between the two functionals is exponentially small, uni- 
formly for r G [y , ^r 1 ] (see (2.12), (2.11) and Proposition 2.2), we get the conclusion 
of Corollary 4. ■ 

Finally, we give in the following some geometric information related to charac- 
teristic points (for the statement in one space dimension, see Theorem 1, Theorem 
2 and the following remark in [22]): 



E (w r (s)) > k(r)E (K ) - C 3 (r )e 



-7s 



(ii) If for some r G [y, -y 1 ] and s > S3, we have 



E (w r (s)) < 2E (k ) - C 3 (r )e 
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Theorem 5 (Geometric considerations on S) 

(i) (Isolatedness of characteristic points) Any r G S R zs isolated. 

(ii) (Corner shape of the blow-up curve at characteristic points) If r G 

S H u>?in k(ro) solitons and Co( r o) G M as center of mass of the solitons' center 
as shown in (4.3) and (4.6), then 

0(r)ue~ 29 ^ o( - r ^ 
nr) + 9(r) ~ - j - ^ee ^ 

| log 1 7™ — 7*0 1 1 2 

r(r)-T(r ) + |r-r | - ■ , ' ^xll, ( 4 - 8 ) 

| log | r — ro 1 1 2 

as r — >• r , where 6{r) = nE^\ an d v = u {p) > 0- 
Proof. See below. 

Remark: As stated in the remark after Theorem 3, our result holds for N = 1, 
with no symmetry assumptions an initial data, for all r G R, even r = 0; when 
N >2, and if G 5, the asymptotic behavior of wq remains open. 
Proof. As in the one-dimensional case with no perturbations 

Remark: Note from (i) that the multi-dimensional version U(x,t) = u(\x\,t) has 
a finite number of concentric spheres of characteristic points in the set {-^ < |x| < R} 
for every R > 1. This is consistent with our conjecture in [22] where we guessed that 
in dimension N > 2, the (N — 1) -dimensional Hausdorff measure of S is bounded in 
compact sets of M. N . Note that this conjecture is related to the result of Velazquez 
who proved in [28] that the (N — l)-dimensional Hausdorff measure of the blow-up 
set for the semilinear heat equation with subcritical power nonlinearity is bounded 
in compact sets of R . 

As a consequence of our analysis, particularly the lower bound on T(r) in (4.8), 
we have the following estimate on the blow-up speed in the backward light cone 
with vertex (r ,T(r )) where r > (for the statement in one space dimension, see 
Corollary 3 in [22]): 

Corollary 6 (Blow-up speed in the backward light cone) For all r > 0, 

there exists 64(70) > such that for all t G [0, T(r )) ; we have 

| log(T(ro) _ () l^ ^ ^ ^ ftWiMTh)-^ 



C4(r )(T(r ) - 1)—^ \x-r \<T(r )-t (T(r ) - 1)^ 

Remark: Note that when ro G 1Z (1 , the blow-up rate of u in the backward light 
cone with vertex (r , T(r )) is given by the solution of the associated ODE u n = u p . 
When r G S PI R+, the blow-up rate is higher and quantified, according to k(r ), 
the number of solitons appearing in the decomposition (4.3). 

Proof. When ro G TZ, the result follows from the fact that the convergence in (3.2) 
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is true also in L°° x L 2 from (3.2) and the Sobolev embedding in one dimension. 
When r G S, see the proof of Corollary 3 of [22] given in Section 3.3 of that paper. 



Proof of Theorems 3 and 5: The proof follows the pattern of the original proof, 
given in [20], [24], [22] and [7]. In the following, we recall its different parts. 

Part 1: Proof of (4.3) without (4.4) nor (4.5) and with the estimate 

Ci+i(s) - Q(s) ->■ oo as s ->> oo (4.9) 

instead of (4.6) (note that (4.9) is meaningful only when k(ro) > 2). 

The original statement of this part is given in Theorem 2 (B) page 47 in [20] and 
the proof in section 3.2 page 66 in that paper. Note that this part doesn't exclude 
the possibility of having k(ro) = or fc(ro) = 1. The adaptation is straightforward. 
As in the non-characteristic case above, one has to use the Duhamel formulation in 
the radial which may be derived from [26]. 

Part 2: Assuming that (4.4) is true, we prove (4.5) with the estimates 

Ms)-Ci(s)\ < C, (4.10) 

T(r)-T(r ) + \r-r Q \ < — — 

| log | r — ro 1 1 2 

instead of (4.6) and (4.8). 

The original statement is given in Propositions 3.1 and 3.13 in [24]. The reader 
has to read Section 3 and Appendices B and C in that paper. The adaptation is 
straightforward, except for the effect of the new terms in equation (2.2), which pro- 
duce exponentially small terms in many parts of the proof (see (3.5)). In particular, 
Lemma 3.11 of [24] has to be changed by adding Ce~ JS where 7 is defined in (2.10) 
to the right of all the differential inequalities. 

Part 3: Proof of (4.4) and the fact that the interior of S is empty. 

The original statement is given in Proposition 4.1 of [24]. The adaptation is as 
delicate as in [25]. In particular, it involves the ruling-out of the occurrence of the 
case where, locally near the origin, the blow-up set of the multi-dimensional version 
U(x,t) is a forward light cone with vertex (0,T(0)) (see Lemma 4.5 page 367 in 
[25]). As in that paper, the proof of the non-occurrence of this case is based in 
particular on a local energy estimate by Shatah and Struwe [26]. For the reader's 
convenience, we adapt in Appendix A that energy estimate to our case (1.2), namely 
to perturbations of the pure power equation (1.9). For the other arguments, we refer 
to the corresponding part in [25] (see Part 3 page 336 in that paper). 
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Part 4: Proof of Theorem 5 with (4.7) and (4.8) replaced by 

— v / r— m — 



„ M , » | (fc(r )-l)(p-l) — ^ i — r | — , s i (fc(r )-l)(p-l) > 

C |log(r-r )| 2 |log(r-r )| 2 

J y ^"Q J I ^* ^"*0 I 

— — - M (fc(r )-i)(P-D ' ^ T(r) - T(r ) + |r - r | <- - - ^pggrip 

C 1 log(r - r ) I 2 I log(r - r ) | 
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The analogous statement in one space dimension with no perturbations is given 
in Theorems 1 and 2 in [22]. Thus, one needs to say how to adapt the analysis of 
the paper [22] to the present case. As in [25], three ingredients are needed in the 
proof: 

- the trapping result stated in Theorem 2; 

- the energy criterion stated in Corollary 4; 

- the dynamics of equation (2.2) around a decoupled sum of solitons performed in 
[22] and presented in Part 3 above. Note that we have already adapted all these 
ingredients to the present context. With this fact, the adaptation given in [25] works 
here. See Part 4 page 371 in that paper for more details. 

Part 5: Proof of (4.6), (4.7) and (4.8) 

This part corresponds to the contributions brought in [7] in the one-dimensional 
case. The orginal statements in the one-dimensional case are given in Theorem 1.1 
and Corollary 1.4 in that paper. Following Part 2 where we proved that (4.3) holds 
with (4.6) replaced by (4.10), a crucial step in one-space dimension was to prove 
that the solitons' centers satisfy the following ODE system for s large enough: 



—Q = e -^r«*-C<-i) _ e -j£i(Ci+i-C0 + q I _!_ ] 



for some r\ > 0. In [7], we were able to use some ODE tools (particularly the 
Lyapunov convergence theorem) to further refine estimate (4.10) and prove that 

Ci(«) = (i( s ) + Co + f^J as s -»■ 00. 

Since we have for all \di\ < 1 and [c^l < 1 

— ^(^2)11^ < C| argtanhrfi — argtanhe^l 

(see estimate (174) page 101 in [20] for a proof of this fact), estimate (4.3) remains 
unchanged if one slightly modifies Q(s) by setting d(s) = Q(s) + Co which is the 
desired estimate in (4.6). That was the argument in one space dimension. 
In our setting, since our perturbative terms contribute with additional exponentially 
decaying terms to the equation (see (3.5) and Part 2 above), we obtain that d(s) 
satisfy the same ODE system (4.11). Thus, the refinements of [7] hold here with no 
need for any further adaptations, and (4.6) holds. 
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As for estimates (4.7) and (4.8), let us point out that in one space dimension, 
they are derived in [7] as direct consequences of (4.6) on the one hand, and on the 
other hand a small improvement of the last argument of the paper [22] based on the 
equation in similarity variables. Since in our setting, (4.6) holds and the equation in 
similarity variables differs from the one dimensional case with exponentially decaying 
terms (see (3.5)), the same argument holds. See Section 2 in [7]. 



4.2 Existence and non-existence of characteristic points 

Proceeding as in [7], we have the following result: 

Theorem 7 (Existence of a solution with prescribed blow-up behavior at 
a characteristic point) For any r > and k > 2, there exists a blow-up solution 
u(r,t) to equation (1.6) with r G S such that 



d s w ro {s) 



( k \ 

^(-l^Md^)) 



i=l 



/ 



— > as s — > oo, 



H 



with 



d,i(s) = -tanhO(s), GO) = GO) + Co 
for some Co £ where Gp) is defined in (4.1). 



(4.12) 



(4.13) 



Remark: When N = 1, we can take ro = 0. When N > 2, the multi-dimensional 
version U(x,t) = u(\x\,t) has a sphere of characteristic points. Note also that this 
result uses the same argument as for Theorem 3, in particular, the analysis of the 
ODE system (4.11). If we simply want an argument for the existence of a blow-up 
solution with a characteristic point without caring about the number of solitons, 
then we have a more elementary proof which holds, however, only when g does not 
depend on \x\. See the remark following Theorem 9 below. 

Remark: Note from (4.13) and (4.1) that the barycenter of GO) * s fixed, in the 
sense that 



GP) + --- + qp) = CiP) + --- + ao) 
k k 



+ Co = Co, Vs> -logT(O). (4.14) 



Note that unlike in the one-dimensional case with a pure power nonlinearity treated 
in [7], we are unable to prescribe the barycenter. Indeed, our equation (1.6) is not 
invariant under the Lorentz transform. 

Remark: We are unable to say whether this solution has other characteristic points 
or not. In particular, we have been unable to find a solution with S exactly equal 
to {0}. Nevertheless, let us remark that from the finite speed of propagation, we 
can prescribe more characteristic points, as follows: 
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Corollary 8 (Prescribing more characteristic points) Let J = {l,...,no} or 
J = N and for all n G J, r n > 0, T n > and k n > 2 such that 



r n J rT n < r n+ i — T n+ i. 



(4.15) 



Then, there exists a blow-up solution u(r,t) of equation (1.2) with {x n \ n G J} C S, 
T{r n ) = T n and for all n G I, 



w Xn {s) 
9 s w Xn (s) 



1 f^(-iy +1 K(d hn ( s )) ^ 



1=1 



/ 



— y as s —7- oo, 



Vz = l,...,k n , di,„(s) = -tanhCi,n(s), 0,n.(s) = + Co,« 
/or some £o,n G R, where Q{s) is defined in (4.1). 

Remark: Again, we are unable to construct a solution with S = {r n | n G J}. 
When TV = 1, we may take r G R. 

Proof of Theorem 7 and Corollary 8: First, note that thanks to condition (4.15) 
which asserts that the sections at t = of the backward light cones with vertices 
(r n ,T n ) do not overlap, Corollary 8 follows from Theorem 7 by the finite speed of 
propagation. As for the proof of Theorem 7, we claim that it follows like in [7], 
since the ingredients of that paper are available here, thanks to the adaptations we 
performed in the previous sections: 

- the analysis of the ODE system (4.11): let us emphasize the fact that we still 
encounter this sytem in our setting. Indeed, that system appears as a projection on 
the null modes of the linearization of equation (2.2) around the sum of decoupled 
solitons, and, as we said in Part 2 page 16, that equation differs from the pure 
power case, only with exponentially small terms (see (3.5)), which are absorbed in 
the O(^) in (4.11); 

- a reduction to a finite dimensional problem: this is done thanks to the 
analysis of the dynamics of the equation in similarity variables (2.2) around the 
sum of decoupled solitons, which we did already for the proof of the isolatedness 
of characteristic points (see (i) of Theorem 5; see [22] for the analysis in one space 
dimension) ; 

- a topolgical argument to solve the finite dimensional problem: this argu- 
ment is based on a different formulation of Brouwer's Theorem. It is independent 
of the equation. 

Note however that one argument of [7] does not work here: the argument that allows 
us to prescribe the barycenter of 0( s )- Indeed, that argument uses the invariance of 
the pure power wave equation under the Lorentz transform, which is no longer the 
case for equation (1.6). ■ 



Let us give in the following a criterion about the existence of characteristic points: 
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Theorem 9 (Existence and generic stability of characteristic points) 
(i) (Existence) Let < a\ < a 2 be two non- characteristic points such that 

w a,i(s) — > 6(ai)K(d ai , ■) as s — > 00 with 6>(ai)6>(a 2 ) = — 1 

for some d ai in (—1,1), in the sense (3.2). Then, there exists a characteristic point 
c G (ai,a 2 ). 

(^(Stability) There exists e > such that if\\(U , Ui) — (U , Ui)\\ R i xL 2 qrjv) < 

6q, then, u(r,t) the solution of equation (1.6) with initial data (uo,ui)(r) = (Uo,Ui)(x) 
if r = \x\ blows up and has a characteristic point c G [01,02]. 

Remark: This statement (valid for N > 2) is different from the original one (The- 
orem 2 in [24]) by two natural small facts: we take positive points a\ and a 2 in (i), 
and we use the multi-dimensional norm in (ii) (of course, from the finite speed of 
propagation, it is enough to take a localized norm instead). When iV = 1, we don't 
need the restriction a\ > 0. 

Remark: If one needs a quick argument for the existence of a blow-up solution for 
equation (1.6) with a characteristic point, then this theorem allows us to avoid the 
heavy machinery of [22], namely the linearization of equation (2.2) around the sum 
of decoupled solitons. Indeed, we have a more elementary argument, based on the 
knowledge of the blow-up behavior at a non-characteristic point on the one hand, 
and on (i) of this theorem on the one hand. However, such an argument uses the fact 
that solutions of the ODE (4.16) associated to (1.6) are also solution to (1.6) and 
this is possible only if g defined in (1.2) does not depend on \x\. For the statement 
with no perturbations, see Proposition 3 page 362 in [25]. For a further justification, 
see the Proof of Theorem 9 below. 

Proof of Theorem 9: As in [25], there is no difficulty in adapting to the present 
context the proof of Theorem 2 of [24] given in Section 2 of that paper, except may 
be for some natural extensions to the radial case. Concerning the second remark 
following Theorem 9, the only delicate point is to find initial data (u , ui) satisfying 
the hypothesis of (i) in Theorem 9. If g does not depend on \x\, then, any solution 
of the ODE 

U'^llir'U + fW+g^O^U'), (4.16) 

is also a solution of the PDE (1.6), and it is enough to take initial data (uq, u{) with 
large plateaus of opposite signs. If g does depend on then this simple idea breaks 
down, and the existence of initial data with characteristic points holds thanks to 
Theorem 7. ■ 

We also have the following result which relates the existence of characteristic 
points to the sign-change of the solution: 

Theorem 10 (Non-existence of characteristic points if the sign is con- 
stant) Consider u(r,t) a blow-up solution of (1.6) such that u(r,t) > for all 
r G (a ,bo) and t < t < T{r) for some real < a < b and t > 0. Then, 
(a ,b ) C K. 
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Remark: When N — 1, we don't need the restriction ao > 0. 
Proof. This result follows from Theorem 3 above exactly as in one space dimension 
with no perturbations (i.e. for equation (1.6) with (/, g) = (0,0)). See the proof of 
Theorem 4 given in Section 4.2 in [24]. ■ 



A A local energy estimate for perturbations of 
the semilinear wave equation 

Let us consider the following perturbation of equation (1.9): 

d 2 t U = AU+ \U\ p - l U + X^f(\^U} + X^g(x\x\, Xt, X~^VU ■ ^- X~^d t U^j , 

(A.l) 

where A > and / and g satisfy (Hf) and (H g ). The equation (A.l) is derived from 
equation (1.2) through the dilation 

A ^ U x (x,t) = Xp= i U(Xx, Xt). 
Using the technique of Shatah and Struwe [26] and introducing 



£(U(t)) 



x <1— t 



dx, (A.2) 



where F is defined in (2.4). We obtain the following local energy estimate 

Lemma A.l (A local energy estimate for perturbations of equation (1.9)) 
For all t G [0, 1), we have 

8{U{t)) < C£(U{0))+C [ [ \U(a,s)\ p+1 dads + CX [ [ \U(x, s)\ p+1 dxds + CX^. 

Jo Jb s Jo J\x\<i-s 



where the lateral boundary is 

B t0 = {(x,t) | < t < t , M = 
Proof. Classical calculation implies 



dt em)) 



' v ' S ■ «, ' 

Li(s) L 2 (s) 

2p f / p+1 X P+l \ 

X~^ / g[\\x\,\t,\~*=tVU • — ,X~^d t U)d t Udx. 
J\x\<i-t ^ fI ' 
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Since |ra-VE/| < |Vf/|, we can say that the term L 2 (s) is negative. By combining 
the estimate |-F(:r)| < C(\x\ + the assumption (H g ) and the fact that the 

terms L±(s) and L 2 (s) are negative, we conclude that, for all A G (0, 1] 

^-S(U(t)) < cf (\fr\U\ + \U\ p+1 )da + CX [ (\£*\dtU\ + \VU\ 2 + \d t U\ 2 )dx. 
at J Bt V / J\x\<i-t v J 

j t E i u ^)) < CX£(U{t))+cJ (A^|C/| + |C/| p+1 )dc7 

+CX I (x^\d t U\ + X^\U\ + \U\ p+l \dx 

J\x\<l~t ^ ' 



So 

d 

Then 



£(U(t)) < CX£(U(t)) +C I \U\ p+1 da + CX [ \U\ p+1 dx + CX^ 

JB t J\x\<l-t 



£{U{t)) < £(U(0))e cxt + C f e cx{t ~ s) [ \U\ p+1 dads 

Jo J B s 

+CX f ' e cxit ~ s) [ \U\ p+1 dxds + CX^ f 

JO J\x\<l-s Jo 



SO 



£{U{t)) < £(U(0))e cxt + C f e cx{t - s) [ \U\ p+1 dads 

Jo J B a 

+CX [ e cx{t - s) [ \U\ p+1 dxds + CX^e cxt 

Jo J\x\<l-s 
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